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ROP-type GTPases of plants function as molecular switches within elementary signal 
transduction pathways such as the regulation of ROS synthesis via activation of NADPH 
oxidases (RBOH-respiratory burst oxidase homolog in plants). Previously, we reported that 
silencing of the Medicago truncatula GTPase MtROP9 led to reduced ROS production 
and suppressed induction of ROS-related enzymes in transgenic roots (MtROP9i) 
infected with pathogenic {Aphanomyces euteiches) and symbiotic microorganisms 
(Glomus Intraradices, Sinorhizoblum meliloti). While fungal infections were enhanced, 
S. meliloti infection was drastically impaired. In this study, we investigate the temporal 
proteome response of M. truncatula MtROP9i transgenic roots during the same microbial 
interactions under conditions of deprived potential to synthesize ROS. In comparison 
with control roots (Mtvector), we present a comprehensive proteomic analysis using 
sensitive MS protein identification. For four early infection time-points (1, 3, 5, 24 hpi), 
733 spots were found to be different in abundance: 213 spots comprising 984 proteins 
(607 unique) were identified after S. meliloti infection, 230 spots comprising 796 proteins 
(580 unique) after G. intraradices infection, and 290 spots comprising 1240 proteins (828 
unique) after A. euteiches infection. Data evaluation by GelMap in combination with a 
heatmap tool allowed recognition of key proteome changes during microbial interactions 
under conditions of hampered ROS synthesis. Overall, the number of induced proteins in 
MtROP9i was low as compared with controls, indicating a dual function of ROS in defense 
signaling as well as alternative response patterns activated during microbial infection. 
Qualitative analysis of induced proteins showed that enzymes linked to ROS production 
and scavenging were highly induced in control roots, while in MtROP9i the majority of 
proteins were involved in alternative defense pathways such as cell wall and protein 
degradation. 
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INTRODUCTION 

RAC/ROP (Rho of plants) are plant-specific small GTPases that 
function as simple binary molecular switches within elementary 
signal transduction pathways by cycling between GTP-bound on 
modes and GDP-bound off modes. In the GTP-bound forms, 
they interact with specific downstream effectors, mediating a wide 
repertoire of molecular stimuli that provoke cellular responses 
(Poraty-Gavra et al., 2013). ROPs integrate many upstream sig- 
nals via the guanine nucleotide exchange factors (GEFs), guanine 
nucleotide dissociation inhibitors (GDIs) and GTPase-activating 
proteins (GAPs), regulating downstream effectors such as Rop- 
interactive CRIB motif-containing proteins (RICs) and interactor 
of constitutively active ROPs (ICRs) (Nagawa et al., 2010). Small 
GTPases are well studied in mammals and yeast cells and have 
been grouped into various subfamilies depending on their func- 
tional properties (Schiene et al., 2000). ROPs are known to 
function in different developmental processes including polarized 



cell growth, pollen tube and root hair development, hormonal 
signaling as well as cell morphogenesis (Yang and Fu, 2007; Liu 
et al, 2010). They are also implicated in regulating several cellular 
processes including vesicle trafficking, cytoskeleton organization 
and dynamics, auxin transport and response to pathogens (Nibau 
et al, 2006; Yalovsky et al, 2008; Lorek et al, 2010; Wu et al, 201 1; 
Poraty-Gavra et al., 2013). The ICRs have also been shown to 
regulate polarized secretion and polar transport of auxin during 
auxin-regulated development. ROP proteins form key regula- 
tory elements for reactive oxygen species (ROS) generation in 
plant cells especially at the plasma membrane by activating the 
NADH oxidases termed as RBOH (for respiratory burst oxidase 
homolog). In Medicago truncatula, 16 putative ROPs were sug- 
gested from assembled EST sequences (Yuksel and Memon, 2008) 
but only seven ROPs have been confirmed (Liu et al., 2010) 
with the rest either being artifacts or redundant. In Arabidopsis 
thaliana, 11 ROPs are identified (Winge et al., 1997). 
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During the initial microbial invasion, the host plant con- 
fers a general defense reaction characterized by rapid activation 
of a wide repertoire of symbiotic or pathogenic defense cellu- 
lar responses. More often, the ROS, especially H2O2, form the 
hallmark of these very early host defense systems that cause 
a hypersensitive reaction culminating in host cell death at the 
site of infection (Puppo et al., 2013). ROS are also said to 
diffuse across the cell membranes via the aquaporins and func- 
tion as second messenger during signal transduction pathway 
hence acting as elementary signal molecules for activation of 
plant defense responses (Borisova et al., 2012). The superox- 
ide and hydrogen peroxide are the typical ROS accumulated in 
host cells. The delicate balance between production and scav- 
enging activity allows the duality in function of ROS to exist 
in plant system orchestrated by a large network of enzymes and 
antioxidative compounds. The scavenging activity of ROS in the 
mitochondria is controlled by the alternative oxidase (AOX), non- 
proton-pumping, alternative type II and the Ca 2+ -dependent 
NADPH dehydrogenase (Steffens et al., 2013). The antioxida- 
tive activity in the cell is provided by molecules such as glu- 
tathione, tocopherols, tannins, phenolic compounds and ROS 
scavenging enzymes such as superoxide dismutase (SOD), ascor- 
bate peroxidase, glutathione peroxidase, catalase as well as other 
non-enzymatic proteins such as metallothioneins and thiore- 
doxin that lead to ROS homeostasis (Steffens et al., 2013). The 
ROS scavenging mechanisms lead to cell wall structural reinforce- 
ment by cross-linking of various extracellular proteins including 
proline-rich glycoproteins to the polysaccharide matrix (Djebali 
et al., 2011). In rice (Oryza sativa), the GTPase OsRacl was 
shown to positively regulate disease resistance by stimulating the 
NADPH-mediated ROS production via direct binding to the cat- 
alytic subunit of NADPH oxidase N-terminal extension specific 
for RBOH proteins (Kawasaki et al, 2006; Jones et al, 2007; 
Nakashima et al., 2008). In vivo fluorescence resonance energy 
transfer (FRET) analysis showed that the Ca 2+ concentration in 
the cytosol may regulate the RBOH-Rac interaction, hence mod- 
ulating the activity of NADPH oxidase in ROS production (Wong 
et al, 2007). 

Legumes (Fabaceae) interact with soil-borne microbes 
(Colditz and Braun, 2010) and are unique in establishing sym- 
biosis with rhizobia bacteria, which ultimately leads to nitrogen 
fixation in the formed structures known as nodules. Through 
their association also with arbuscular mycorrhizal fungi, legumes 
benefit by acquiring macronutrients including phosphorus and 
nitrogen, as well as most likely an array of micronutrients in 
exchange for up to 20% of the plant-fixed carbon (Finlay, 2008). 
ROPs are shown to play a key role during rhizobia infection in 
the process of nodule development and also during the establish- 
ment of mycorrhizal association (Berken, 2006). Expression of 
MtROP3, MtROP5 and MtROP6 in M. truncatula increased after 
rhizobia inoculation as reported by Liu et al. (2010), while in 
L. japonicus, the LjROP6 was shown to act as a positive regulator 
of infection thread formation during rhizobia infection (Ke et al., 
2012). 

Previously, we reported that silencing of MtROP9 impairs 
rhizobial infection but positively regulates root colonization 
by arbuscular mycorrizal fungi G. intraradices and oomycete 



pathogen A. euteiches (Kiirika et al., 2012). The infection process 
in MtROP9i transgenic roots was characterized by clear reduction 
of ROS accumulation in the cells, and by marked transcriptional 
suppression of ROS-related enzymes such as RBOH and catalase. 
Both symbiotic and pathogenic interactions are known to induce 
oxidative burst coupled with induction of defense-related prod- 
ucts such as PR proteins, where the difference in the two forms 
of interactions is suggested to be of quantitative in nature espe- 
cially with regard to ROS production. In addition, the generation 
of ROS suppresses the expression of PR genes (Peleg-Grossman 
et al., 2012). In M. truncatula, the MtSpkl gene encoding a 
putative protein kinase was induced by exogenous application 
of H2O2 as well as nodulation factor indicating the functional 
role of ROS in regulating genes directly linked to rhizobia sym- 
biosis (Andrio et al., 2013). Transient decrease in MtRBOHs 
gene expression was reported to lead to decrease in ROS efflux 
observed 1 h after M. truncatula roots treated with NF (Lohar 
et al, 2007). 

In this study, we have utilized the previously investigated 
sequence of MsRacl sequence from Medicago sativa for RNA 
interference (RNAi) -mediated gene silencing in the model legume 
M. truncatula where we identified a M. truncatula sequence 
ortholog annotated as MtROP9 (TC173331; Dana-Farber Cancer 
Institute M. truncatula Gene Index [MtGI]; Quackenbush et al., 
2001). A gene-specific region of MtROP9 was selected for RNAi 
gene knockdown with Agrobacterium rhizogenes used as a vector. 
Evaluation of MtROP9i root proteome maps via 2D IEF SDS- 
PAGE and MS after symbiotic and pathogenic interactions at the 
very early timepoints of infections revealed changes in protein 
profiles as clearly visualized using the heatmap-GelMap tool. 

MATERIALS AND METHODS 
CONSTRUCTION OF RNAi VECTOR 

Transgenic MtROP9i roots were produced as reported 
previously (Kiirika et al., 2012), using the binary vector 
pK7GWIWG2(II)::DsRED (kindly provided by E. Limpens; 
Limpens et al., 2005) containing the gene for red fluorescent 
marker DsREDl. The vector was modified by insertion of two 
sequence cassettes (in the sense-antisense direction) encoding 
parts of the putative effector (G2) and GTPase (G3) domains of 
the MsRacl ortholog MtROP9 (TC173331; GenBank accession 
no. AF498359). Binary vectors for gene knockdown by RNAi 
were constructed using the Gateway technology (Invitrogen Life 
Technologies). Gene-specific oligonucleotides (ROP9attbl_for, 
5'-attBl GTGTTACTGTTGGTGATG-3';ROP9attb2_rev,5'-attB 
2-ACGCCTTCACGTTCTCC-3') with attached attB adapters 
were obtained from the Medicago sativa MsRacl sequence 
(GenBank accession no. AJ251210; Schiene et al, 2000). Using 
these oligonucleotides, amplification of a 461 -bp fragment from 
Medicago truncatula cDNA was carried out and the PCR products 
purified using the QIAquick PCR purification kit (Qiagen) 
and cloned into the pDONR221 donor vector (Invitrogen). In 
the second cloning step, the inserts were transferred into the 
Gateway-compatible binary vectors mentioned above followed 
by transformation into Agrobacterium rhizogenes strain Arqual 
(Quandt et al, 1993) using standard methods. The presence 
of the MtROP9 sequence fragments and the sense-antisense 
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orientation of the cloned fragments in the T-DNA were con- 
firmed by sequencing the construct, while non-modified binary 
vectors were transformed into A. rhizogenes Arqual as a control. 

GENERATION OF TRANSGENIC ROOTS AND INFECTION ASSAYS 

M. truncatula MtROP9i and Mtvector transgenic roots were 
generated according to the Boisson-Dernier et al. (2001) using 
M. truncatula (Jemalong A17) wild-type plantlets as described 
before (Colditz et al, 2007). M. truncatula composite plants 
with roots transformed by A. rhizogenes were cultured stably 
onM medium (Becard and Fortin, 1988) containing 25 mg L21 
kanamycin for selection and decreasing concentrations of 350 to 
Omg L21 ticarcillin disodium/clavulanate potassium (Duchefa) 
to stop growth of A. rhizogenes. Twelve individual populations 
of composite plants with MtROP9i and Mtvector transgenic 
roots were generated independently via A. rhizogenes transfor- 
mation, containing at least 200 plants each. Populations were 
cultured onM medium and kept in the growth chambers at 
22°C, 65% humidity, 16-h photoperiod at 220 |xE mT 2 s _1 . 
Inoculation with A. euteiches (ATCC 201684) was carried out as 
described before (Colditz et al., 2007). Each transgenic root pop- 
ulation was infected with 500 mL of lake water containing 250,000 
A. euteiches vital zoospores. Inoculation with G. intraradices was 
performed using commercially available inoculums (Granular 
AMF inoculum; BIORIZE). Infection with Sinorhizobium meliloti 
wild-type strain Rm2011 was performed as described previously 
(Schenkluhn etal., 2010). 

PROTEIN ISOLATION, 2D IEF SDS-PAGE AND GEL EVALUATION 

Phenol extraction of total protein from the cells was carried out 
as described previously (Colditz et al, 2004, 2005, 2007). For IEF, 
3 mg of protein was diluted with 350 |xl of rehydration buffer, 
consisting of 8 M urea, 2% (wt/vol) CHAPS, 100 mM dithiothre- 
itol, 0.5% (vol/vol) IPG buffer for correspondent pi range (pH 
3-11 non-linear [NL]; Amersham Pharmacia Biotech, Uppsala, 
Sweden), and a trace of bromphenol blue. 2D IEF SDS-PAGE 
was performed for the three infections and at four timepoints of 
harvesting by combining the IEF strips (IPGphor system) with a 
sodium dodecyl sulfate-tricine gel electrophoresis (Protean II XL, 
20 by 20 cm; BioRad, Richmond, CA, U.S.A.) as already described 
(Colditz et al, 2005). Gels were stained with 0.1% (wt/vol) 
Coomassie Brilliant Blue (BioRad) overnight and scanned on a 
UMAX Power Look III Scanner (UMAX Technologies, Fremont, 
CA, U.S.A.). Gels were evaluated using Delta 2D software, ver- 
sion 4.0 (Decodon, Greifswald, Germany) with three replicates 
per group (1, 3, 5, and 24 hpi). Spots detection was done automat- 
ically and occasionally corrected manually. In gel normalization 
was performed using the Delta 2D software for the overlays of 
three replicate gels each. Spots with a relative spot volume of less 
than 0.05% were deleted and the significant abundance of spots 
between MtROP9i and Mtvector groups was determined using 
Student's f-test (confidence interval >95%) based on the relative 
spot volume. 

VISUAL EVALUATION OF PROTEIN INDUCTION PATTERNS VIA GELMAP 
EXTENDED BY A HEATMAP TOOL 

To explicitly present the voluminous protein dataset comprising 
of all significantly induced proteins from the three infections, i.e., 



S. meliloti, G. intraradices and A. euteiches in both transgenic root 
populations MtROP9i and Mtvector, the novel software GelMap 
(http://www.gelmap.de) which is utilized for protein annotation 
was used (Klodmann et al., 2011; Senkler and Braun, 2012). In 
addition, a heatmap tool was integrated to the GelMap mod- 
ule allowing a clear visualization of induction patterns of all 
identified proteins according to differences in abundance. This 
technique utilizes the inbuilt function filters in the GelMap 
software for any user-defined cluster of choice, based on spe- 
cific physiological functions of proteins. The presented Heatmap 
(Figure 2; https://gelmap.de/532) was created by clustering all 
protein induction values of proteins of similar physiological func- 
tions both for major and sub-categories found at a certain selected 
infection time-point. The total inductions per sub-category pre- 
sented as individual values were automatically filtered and dis- 
played on a matrix. The Heatmap was then generated by applying 
a color-coding system which visually indicates the range of pro- 
tein induction for a certain (sub-) category from dark-red to 
light-red, corresponding to high protein inductions and low pro- 
tein inductions, respectively. Proteins with the highest levels of 
induction either in MtROP9i or Mtvector were selected or ranked 
first as the most predominant physiological categories in the 
Heatmap-GelMap (Table 2). 

MASS SPECTROMETRY AND CREATION OF HEATMAP-GELMAP 

Protein spots of 1.4 mm diameter were cut from Coomassie 
stained gels using a GelPal Protein Excision manual spot picker 
(Genetix, Great Britain) and in-gel digested with Trypsin as 
described by Klodmann et al. (2011). Tryptic peptides were fur- 
ther analyzed by nanoHPLC (Proxeon, Thermo Scientific) cou- 
pled to electrospray ionization quadrupole time of flight MS 
(micrOQTOF Q II, Bruker Daltonics), using all settings and 
parameters as described previously (Klodmann et al., 2011). 
Data processing and protein identification was carried out with 
ProteinScape 2.0 (Bruker Daltonics) and the MASCOT search 
engine querying three Medicago -specific protein databases [Mf3.5 
ProteinSeq, NCBI Medicago truncatula protein, and Mtf(asta)2] 
available at the LegProt db (Lei et al., 2011) as well as SwissProt, 
using the following parameters: trypsin/P; one missed cleavage 
allowed; fixed modifications: carbamidomethylation (C), vari- 
able modifications: acetylation (N) and oxidation (M); precursor 
ion mass tolerance, 30ppm; peptide score >24; charges 1C, 
2C, 3C. Protein and peptide assessments with MASCOT scores 
above 25 were considered. Heatmaps were created using the 
total significant induction (>1.5 fold) of proteins per physiolog- 
ical function as annotated via Swissprot. Induction of proteins 
with similar functions was clustered and presented as values, 
which were filtered and displayed on a matrix. Color-coding 
scheme ranging from dark red (high induction) to light red (low 
induction) was applied. A Heatmap image (.jpg) was assigned 
x- and y-coordinates to specifically allocate total inductions of 
each protein category into a corresponding color-code range on 
a Heatmap matrix (Figure 2). The Heatmap was linked with 
the reference GelMap platform (Senkler and Braun, 2012). An 
Excel (Microsoft) file containing all protein information includ- 
ing MS/MS results and the corresponding Heatmap image (.jpg) 
were then imported into GelMap portal. Information on GelMap 
creation can be accessed at http://www.gelmap.de/howto. The 
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heatmap tool integrated to the GelMap, as utilized for protein 
annotation (Klodmann et al., 2011), allowed explicit presenta- 
tion of voluminous dataset and clear visualization of induction 
patterns of all identified proteins according to differences in abun- 
dance. This technique utilizes the inbuilt function filters in the 
GelMap software for any user-defined cluster of choice, based on 
specific physiological functions of proteins. It clusters all protein 
induction values of proteins of similar physiological functions 
both for major and sub-categories found at a certain selected 
infection time-point. 

RESULTS AND DISCUSSION 
EXPERIMENTAL SETUP 

The aim of this study was to evaluate the temporal proteome 
response of M. truncatula roots to symbiotic and pathogenic 
infection under conditions of significantly suppressed poten- 
tial for the legume host to synthesize ROS as a common and 
early defense mechanism. We generated transgenic M. truncat- 
ula (Jemalong A17) root populations expressing an RNAi hairpin 
construct for the RAC-type GTPase MtROP9 (TC173331; MtGI 
at Dana-Farber Cancer Institute; Quackenbush et al., 2001) as 
described and characterized by us previously (Kiirika et al., 2012). 
The transgenic roots termed MtROP9i revealed significantly 
reduced levels of ROS production as compared to transgenic 
empty vector control (Mtvector) roots (Kiirika et al., 2012). As 
a consequence, MtROP9i roots exhibited clearly altered infection 
levels when inoculated with symbiotic and pathogenic microbes 
(Kiirika et al, 2012). 

Twelve individual populations of MtROP9i were constructed 
independently via A. rhizogenes-mediated root transforma- 
tion, each representing at least 200 stably growing composite 
plants with transgenic roots. For our infection assays and pro- 
teomic analyses, only MtROP9i and Mtvector roots compris- 
ing on average 60% and more of transformed (transformation 
marker/DsRED-positive) roots were considered. For validation of 
the effective knockdown in MtROP9i transgenic roots, MtROP9 
transcript abundance was determined by reverse transcription 
(RT)-quantitative PCR as previously described (Kiirika et al., 
2012). The relative expression of MtROP9 was drastically reduced 
about >90% in these root populations and did not increase 
after microbial infection when compared with Mtvector roots 
(data not shown). Importantly, MtROP9i transgenic roots of all 
considered populations revealed no significant ROS production 
nor its accumulation after infections with the chosen compati- 
ble root microbes as confirmed via in planta luminometric ROS 
assays while Mtvector roots did (data not shown). Inoculation 
assays were conducted using rhizobial bacterium S. meliloti as well 
as arbuscular mycorrhizal fungus G. intraradices for symbiotic 
interactions, and the legume root pathogen A. euteiches (oomy- 
cota). The microbial inoculations resulted in similar infection 
patterns for MtROP9i and Mtvector transgenic roots as previ- 
ously described (Kiirika et al., 2012) and thus are not explicitly 
shown here. For the protein analyses, proteins from four inde- 
pendent isolations and from four different time-points of harvest 
(1, 3, 5, and 24hpi) were separated using two-dimensional (2-D) 
IEF/SDS PAGE. Three of the proteome maps achieved per sin- 
gle investigation were selected for Delta 2D analysis in order to 



decipher protein spots with significantly differential abundance 
(Figure 1). In parallel, proteins of interest revealing differences in 
abundance were excised from the gels for tryptic digestion and 
MS-based protein identification. 

PR0TE0MIC PROFILING OF MtR0P9i AND Mtvector TRANSGENIC 
ROOTS DURING SYMBIOTIC AND PATHOGENIC INTERACTIONS 

To identify changes in the proteomes of M. truncatula MtROP9i 
transgenic roots infected with S. meliloti, G. intraradices and 
A. euteiches, proteome maps of four early infection points ( 1 h, 
3 h, 5 h, and 24 h) were prepared. Three Coomassie-stained 
gels were produced of infected root tissue of three individual 
MtROP9i and Mtvector populations for each time-point. These 
gels were evaluated by the Delta2D (v4.4) software (Decodon 
GmbH, Greifswald, Germany). Using student's f-test with a confi- 
dence interval >95%, significant changes in spot pattern depicted 
as differences in spot abundance (>1.5 fold) for different infec- 
tions were determined based on the spot volume. These spots 
were selected for further MS analyses. For this purpose, the 
spots were manually excised from the Coomassie-stained gels, 
then in-gel digested with Trypsin and analyzed via nLC ESI-MS. 
Protein identification was carried out based on the Medicago- 
specific protein databases from the LegProt db (Lei et al., 2011), 
which allowed high rates of protein identification. The protein 
search tool ProteinsScape 2.0 (Bruker Daltonics) and MASCOT 
search engine were used for querying the three Medicago-specific 
databases [Mt3.5 ProteinSeq, NCBI Medicago truncatula protein 
and Mtf (asta) 2 ] as well as Swissprot. 

By following this procedure, in total 733 spots were identified 
to be different in abundance at the four considered time-points 
in gels of infected MtROP9i roots as compared to control roots. 
213 spots from S. meliloti infections were found with different 
abundances, representing a total of 984 MS-identified proteins 
that comprised 607 unique proteins. From the 984 proteins, 385 
and 568 proteins were identified in MtROP9i and the control, 
respectively. A comparable number of 230 spots of differences 
in abundance were found in gels for the G. intraradices infec- 
tions, whereof 796 proteins were identified revealing a total of 580 
unique proteins. From these 796 proteins, 311 and 485 proteins 
were induced in MtROP9i and the control, respectively. Infections 
with A. euteiches revealed 290 spots with different abundances 
after the software-based evaluation, giving an increased number 
of 1240 proteins in total, whereof 828 were unique proteins. 1240 
proteins consist of 456 and 784 proteins induced in MtROP9i and 
the control, respectively, after pathogenic infection. 

Overall, the total number of induced proteins identified in 
MtROP9i transgenic roots after microbial infections was found 
to be lower as compared to the protein number in Mtvector. 
This can be attributed to the knock-down of the signaling pro- 
tein MtROP9 which is directly involved in several early signaling 
cascades. However, the highest individual protein induction levels 
were found in Mtvector roots infected with S. meliloti as com- 
pared to fungal infections with A. euteiches and G. intraradices 
(Figure 1; Tables 1,2). 

Considering the induction pattern along a timeline for the 
three infections, protein induction was high at the early points 
of infection in Mtvector with majority of proteins reaching 
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MrROP9i and Mtvector with S. meliloti infection 



protein SGT1 homolog, 159, 3h [17\ 
SOD[Cu-Zn),71, 24h[511 
prohibitin, 246, 24h 
germin-like protein subfamily 2 member 3, 305, 24fy 
[431 

RBOHB, 410, 5h [23] 
germin-like protein subfamily 1 member 1, 457, 5h 
[41]" 

isopentenyl pyrophosphate isomerase, 851, 24h, 
[45] 

small heat shock protein, chloroplastic, 777, 241) 
[78] 

PR10-1 protein, 775, 3h [521 
glutathione S-transferase, 793, 5h [33J 




dihydroflavanol-4-reductase 1, 266, 3h 
[34] 

RBOHB, 411, 3h [33] 

peroxidase, 426, 5h [84] 

RBOHI,464,5h [55] 

C BL-interacting protein kinase 15 , 417, 

3h [56] 

nodule-specific cysteine-rich peptide 96, 
464, 5h [55] 

auxin responsive SAUR protern,718, 5h [47] 
pectinesterase, 738, 5h [62] 
PR10-1 protein, 772, 3h [52] 
glutathione S-transferase, 796, 5h [53] 



MtROP9i and Mtvector with G. intraradices infection 



Pprg2 protein, 10, 5h [61 
classlO PR protein 1, 6, 24h [78J 
haem peroxidase,plant/fungal/bacterial,52,3h[5] 

soluble inorganic pyrophosphatase, 402,24h[12J 

endochitinase, 408, 5h [3] 

inositol-3-phosphate synthase, 422, 5h [23] 

GTP-binding protein-related, 496, 24h [13J 

NAD(P)H-dependent6'-deoxychalcone synthase, 
422,24h[24] 

clathrin assembly protein:At4g32285,707, 3h [4J 



MLP-like protein, 777, 3h [41 




J4-3-3-like protein, 477, 3h [34] 
cationic peroxidase 1, 219, 3h [43] 
y-H(+}-ATPase subunit A, 112c, 5h [39] 

,peroxidase2, 262, 3h [23] 

blue (type 1) copper domain, 545, 3h [12] 



polygalacturonase inhibitor protein, 570, 
5h [23] 

progesterone 5- beta- reductase,725, 

3h [16] 

PR protein 4A, 747, 24h[20] 
PR protein la, 862, 5h [13] 
glutathione peroxidase, 761, 5h [23] 



MtROP9i and Mtvector control with A. euteiches infection 



nectarin IV, 24, 3h [10], 

ABA-responsive protein ABR17, 42, 5h [71 

isoflavonoid xylocosyltransferase, 913, 5h [4J 
rhicadhesin receptor, germin-like protein, 391, 
3h [10] 

inositol-3-phosphate synthase, 422, 5h [4j 

protein P21 kinase inhibitor, 615, 3h [50J 
PR 5 protein, 625, 3h [12J 
glycine-rich RNA-binding, ABA-inducible protein* 
678, 24h [3] 

nascent polypeptide-associated complex beta, 
subunit, 729, 5h [29] 

glutathione S-transferase 796,3h[30] , 



-*k& • - - 



A M P I o 



polygalacturonase inhibitor protein, 34, 
3h [14] 

glutathione S-transferase, 103, 5h [50] 
RBOHI,96,5h[21] 

translationally-controlled tumor protein 
homolog, 358, lh [5] 
germin-like protein subfamily 2 
member 3, 368, 5h[7] 
peroxidase i, 336, 5h [21] 
hydroxycinnamoyl-CoA quinate hydroxy- 
cinnamoyltransferase, 845, 3h [25] 

falreticulin-1, 578, 5h [7] 



PR10-1 protein, 760, 3h [50] 
50D[Cu 2n],779,5h [46] 



FIGURE 1 | Proteins of high abundance in MtROP9i (red spots) and 
Mtvector (green spots) after S. meliloti (A), G. intraradices (B) and 
A. euteiches (C) infections for the predominant physiological categories 
defense response, stress response, signal transduction, secondary 



metabolite biosynthesis and transport. Proteins of equal abundance 
appear as white spots. The protein name, spot number, hours post 
inoculation (h) is given per protein as well as the level of abundance based 
on > 1.5-fold, shown in square brackets. 
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Table 1 | Protein categories induced only in specific infections. 

Physiological function Physiological function (sub-category) 

(main category) 



S1NORHIZOBIUM MELILOTI 



Secondary metabolite 


Aromatic compounds biosynthesis 


biosynthesis 




Signal transduction 


Hormone metabolism (auxin responsive 




SAUR protein) 




Nodulation 


Nucleic acid metabolism 


Pyrimidine biosynthesis 


Energy metabolism 


Sulfur metabolism 




Pyruvate and TCA cycle metabolism 




(malate dehydrogenase) 


APHANOMYCES EUTEICHES 


Defense response 


Chaperone activity 




Protease inhibitor (kunitz-type) 




Proteolysis 


Ion binding & cofactor 


Oxidoreductase activity 


activity 




Lipid/fatty acid 


Lipid & fatty acid biosynthesis 


metabolism 




Protein folding & 


Protein modification 


processing 






RNA binding 


Signal transduction 


Protein translocation 


Stress response. 


Hormone metabolism 


GLOMUS INTRARADICES 


Protein folding and 


Protein modification 


processing 




Signal transduction 


Protein-protein interaction(14-3-3-like 




protein) 




Ion binding(blue (Type 1 ) copper domain) 




Intracellular components binding (Ran 




binding protein) 


Stress response 


Heat shock protein (stromal 70 kDa heat 




shock-related protein) 


Energy metabolism 


Oxidative phosphorylation (ATP synthase 




subunit delta) 


Lipid/fatty acid 


Lipid & fatty acid biosynthesis 


metabolism 




Amino acid metabolism 


NAD binding 


Transport 


Protein transport 



List of physiological function categories (main and sub-categories) of all proteins 
identified exclusively upon specific infections with S. meliloti, A. euteiches and 
G. intraradices. 



maximum induction even at 3 hpi as compared to MtROP9i. 
In MtROP9i, low protein induction at the early infection time- 
points much likely indicates the absence of the first line of defense 
comprising of ROS defense signaling. Thus, early cellular defense 
reactions were affected by the silencing of the signaling protein 
MtROP9. In addition, also the majority of other alternatively 
induced enzymes which are not directly related to ROS signal- 
ing such as those involved in cell wall and protein degradation 
as well as PR proteins were highly induced especially during the 
advanced stages of infection. 



CATEGORIZATION OF INDUCED PROTEINS IN MtR0P9i AND CONTROL 

All identified proteins were further classified according to their 
physiological functions and ordered based on the most predom- 
inant physiological categories, with four categories representing 
the majority of induced proteins: (i) defense response, (ii) stress 
response, (iii) signal transduction, and (iv) secondary metabolite 
biosynthesis. 

In Figure 1, proteins with the highest differential abundance 
in each infection were filtered and visualized on a two-colored 
gel image channel. Evaluation of protein abundances revealed 
suppressed induction of ROS-related enzymes in MtROP9i trans- 
genic roots as compared to Mtvector roots. In this category, 
18 and 43 proteins were differentially induced in MtROP9i and 
Mtvector, respectively, after S. meliloti infection, while 29 and 42 
proteins were induced in MtROP9i and Mtvector, respectively in 
G. intraradices infection. In A. euteiches infection, 31 and 69 pro- 
teins were induced in MtROP9i and Mtvector, respectively. These 
proteins are intuitively linked to ROS production, hypersensitive 
response reaction as well as enzymes responsible for ROS scaveng- 
ing. They include RBOH1 (55 fold at 5 h) and peroxidase (84 fold 
at 5 h) induced after S. meliloti infection in Mtvector as compared 
to MtROP9i roots. SOD [Cu-Zn] (46 fold at 5 h), RBOHl(21 fold 
at 5 h), peroxidase 1 (25 fold at 3 h) and peroxidase pxdc (7 fold 
at 5 h) were induced after A. euteiches infection as well as cationic 
peroxidase 1 (43 fold at 3h), glutathione peroxidase (23 fold at 
5h), peroxidase 2 (23 fold at 3h) and thioredoxin hi (22 fold 
at 5h) after G. intraradices infection in Mtvector compared to 
MtROP9i transgenic roots. 

Proteins that were found be highest in abundance after 
S. meliloti infection include SGT1 homolog (77 fold at 3 h), pro- 
hibitin (73 fold at 24 h), germin-like protein subfamily 2 member 
3 (43-fold at 24 h) and isopentenyl pyrophosphate isomerase 
(45 fold at 24 h) in MtROP9i roots, while pectinesterase (62 
fold at 5h), peroxidase (84 fold at 5h), nodule-specific cysteine- 
rich peptide 96 (78 fold at 3h) and dihydroflavanol-4-reductase 
1 (34 fold at 3h) were found in Mtvector roots (Table 2). In 
G. intraradices infections, PR10-1 (6 fold at 5h), inosito-3- 
phosphate synthase (23 fold at 5 h) and NAD(P)H-dependent 
6'-deoxychalcone synthase (24 fold at 24 h) were identified in 
MtROP9i roots while 14-3-3-like protein (34 fold at 3 h), cationic 
peroxidase 1 (43 fold at 3 h) and progesterone 5-beta- reductase 
(16 fold at 3h) in Mtvector roots. In A. euteiches infection, 
protein P21 kinase inhibitor (50 fold at 3h), isoflavonoid xylo- 
cosyltransferase (4 fold at 5h), and ricadhesin receptor, germin- 
like protein (9 fold at 3 h) were identified in MtROP9i while 
PR10-1 protein (50 fold 3h), peroxidase 1(25 fold at 3h), and 
hydroxycinnamoyl-CoA quinate hydroxy-cinnamoyltransferase 
(21 fold at 5 h) in Mtvector roots (Figure 1; Table 2). 

OVERVIEW OF PROTEIN INDUCTION VIA HEATMAP 

To explicitly present the voluminous protein dataset compris- 
ing of all significantly induced proteins from the three infec- 
tions applied in both transgenic root populations MtROP9i and 
Mtvector, the novel online software GelMap (http://www.gelmap. 
de) which is utilized for protein annotation was used (Klodmann 
et al., 2011; Senkler and Braun, 2012). In addition, a novel 
Heatmap tool was integrated to the GelMap module allowing a 
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Table 2 | List of proteins with highest induction in abundance for selected physiological categories. 



Spot No. 


hpi 


Protein name 


Physiological 
function 


Physiological function 
sub-category 


Diff. ab. 


Accession No. 




MtROP9i_A EUTEICHES 












796 


3R 


Glutathione S-transferase 
GST 


Defense response 


ROS detoxification & 
degradation 


30 


MTR_5g090910 


729 


5R 


Nascent 

polypeptide-associated 
complex subunit beta 


Defense response 


Chaperone activity 


29 


Medtr3g020520.1 


625 


3R 


Pathogenesis-related 
thaumatin-like protein 


Defense response 


Pathogenesis-related 
protein 


12 


Medtr5g01 0640.1 


24 


3R 


Nectarin IV 


Defense response 


Proteolysis 


10 


Q3KU27_NICLS 


75 


5R 


ABA-responsive protein 
ABR17 


Defense response 


Pathogenesis-related 
protein 


7 


ABR17_PEA 


615 


3R 


Protein P21 kinase inhibitor 


Stress response 


Signal transduction 
inhibitor 


50 


P21_SOYBN 


223 


3R 


Cold shock protein-1 


Stress response 


Zinc ion binding 


50 


Q8LPA7_WHEAT 


772 


5R 


Thioredoxin fold 


Stress response 


Antioxidative and 
peroxidase activity 


4 


Q06H32_ARAHY 


678 


24R 


Glycine-rich RNA-binding, 
abscisic acid-inducible 
protein 


Stress response 


RNA binding 


3 


GRPA_MAIZE 


285 


3R 


Peroxidase 1 


Stress response 


Antioxidative and 
peroxidase activity 


3 


Q9XFL3_PHAVU 


755 


5R 


BiP 


Signal transduction 


ATP & nucleotide binding 


16 


Q587K1_SOYBN 


391 


3R 


Rhicadhesin receptor, 
germin-like protein 


Signal transduction 


Carbohydrate 
metabolism 


9 


Medtr7g111240.1 


22 


3R 


Monocopper oxidase-like 
protein SKU5 


Signal transduction 


Ion binding 


5 


Medtr4g101 650.1 


422 


5R 


lnositol-3-phosphate 
synthase 


Signal transduction 


Protein-protein 


4 


IN01_NICPA 








interaction 






271 


24R 


Receptor-like protein 
kinase 


Signal transduction 


Protein-protein 
interaction 


3 


MTR_3g009050 


913 


5R 


Isoflavonoid 
xylocosyltransferase 


Secondary 
metabolite 
biosynthesis 


Phytoalexin biosynthesis 


4 


Q7XZD0_GLYEC 


379 


1R 


VDAC1.3 


Transport 


Ion transport 


3 


TC179231 


Mtvector_A EUTEICHES 












760 


3V 


PR10-1 protein 


Defense response 


Pathogenesis-related 
protein 


50 


MTR_2g035100 


103 


5V 


Glutathione S-transferase 


Defense response 


ROS detoxification & 
degradation 


50 


GSTF5_ARATH 


779 


5V 


SOD[Cu-Zn] 


Defense response 


ROS detoxification & 
degradation 


46 


SODC_PEA 


96 


5v 


RBOH I 


Defense response 


ROS production & 
hypersensitive response 

i nrl i ir i tir\n 

1 1 IU ULi LIU] 1 


21 


RBOHLARATH 


34 


3V 


Polygalacturonase inhibitor 
protein 


Defense response 


Cell wall degradation 


14 


MTR_7g023690 


845 


3V 


Peroxidase 1 


Stress response 


Antioxidative and 
peroxidase activity 


25 


PER1_ARATH 


29 


3V 


Thioredoxin h 


Stress response 


Antioxidative and 
peroxidase activity 


8 


A1BLP6_MEDTR 


305 


5V 


Peroxidase pxdC 


Stress response 


Antioxidative and 
peroxidase activity 


7 


Q40366_MEDSA 


572 


3V 


HSP70-related protein (BiP) 


Stress response 


Chaperone activity 


6 


Q587K1_SOYBN 



(Continued) 
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Table 2 | Continued 

Spot No. hpi Protein name Physiological Physiological function Diff. ab. Accession No. 

function sub-category 



391 


24V 


Peroxidase2 


Stress response 


Antioxidative and 
peroxidase activity 


5 


O24080_MEDSA 


29 


3V 


Putative membrane protein 


Signal transduction 


Protein-protein 


8 


YCF1_LOBMA 




ycf! 




interaction 






578 


5V 


Calreticulin-1 


Signal transduction 


Calcium ion second 
messenger 


7 


A0A762_SOYBN 


368 


5V 


Germin-like protein 
subfamily 2 member 3 


Signal transduction 


Carbohydrate binding 


7 


GL23_ARATH 


556 


3V 


Endoplasmic reticulum 
HSC70-cognate binding 
protein 


Signal transduction 


ATP & nucleotide binding 


6 


022639_SOYBN 


358 


1V 


TCTP 

(translationally-controlled 


Signal transduction 


Calcium ion second 
messenger 


5 


TCTP_MEDSA 




tumor protein homolog) 










336 


5v 


Hydroxycinnamoyl-CoA 
quinate hydroxycinnamoyl- 


Secondary 
metabolite 


Phenylpropanoid 
biosynthesis 


21 


TC 194348 




transferase 


biosynthesis 








391 


24V 


Short-chain dehydrogenase 


Transport 


Intracellular transport 


5 


TIC32_ARATH 




TIC 32, chloroplastic 










MtROP9i. 


S. MELILOTI 












159 


3R 


Protein SGT1 homolog 


Defense response 


Protein degradation 


77 


SGT1_ORYSJ 


71 


24R 


SOD [Cu-Zn] 


Defense response 


ROS detoxification & 
degradation 


51 


MTR_7g114240 


793 


5R 


Glutathione S-transferase 


Defense response 


ROS detoxification & 
degradation 


33 


GSTUM_ARATH 


775 


3R 


PR10-1 protein 


Defense response 


Pathogenesis-related 
protein 


52 


MTR_2g035100 


410 


5R 


RBOH B 


Defense response 


ROS production & 


33 


RBOHB_SOLTU 


246 


24R 


Prohibitin 


Stress response 


hypersensitive response 
induction 

Chaperone activity 


73 


Medtr3g008250.1 


777 


24R 


Small heat shock protein, 
chloroplastic 


Stress response 


Heat shock protein 


78 


HS21C_Wheat 


894 


5R 


Peroxidase pxdC 


Stress response 


Antioxidative and 
peroxidase activity 


72 


Q40366_MEDSA 


408 


5R 


Thioredoxin superfamily 
protein 


Stress response 


Antioxidative and 
peroxidase activity 


23 


A0JQ12_ARATH 


623 


3R 


Heat shock 70 kDa protein, 
mitochondrial 


Stress response 


Chaperone activity 


22 


HSP7M_PEA 


305 


24R 


Germin-like protein 
subfamily 2 member 3 


Signal transduction 


Carbohydrate binding 


43 


GL23_ARATH 


457 


5R 


Germin-like protein 


Signal transduction 


Carbohydrate binding 


41 


GL11_ARATH 


543 


5R 


subfamily 1 member 1 


Signal transduction 


Calcium ion second 
messenger 


29 


TCTP_MEDSA 


TCTP 

(translationally-controlled 




tumor protein homolog) 










851 


24R 


Isopentenyl pyrophosphate 
isomerase 


Secondary 
metabolite 
biosynthesis 


Terpenoid metabolism 


45 


Q9AVG8_TOBAC 


895 


24R 


v-type proton ATPase 
subunit B1 


Transport 


Ion transport 


29 


VATB1_ARATH 




738 


5V 


Pectinesterase 


Defense response 


Cell wall degradation 


62 


MTR_8g 104620 


464 


5V 


RBOH I 


Defense response 


ROS production & 
hypersensitive response 
induction 


55 


RBOHLARATH 



(Continued) 
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Table 2 | Continued 



Spot No. 

796 
775 
410 

426 



hpi Protein name 

5V Glutathione S-transferase 

3V PR10-1 protein 

3V RBOH B 

5V Peroxidase 




Physiological 
function 

Defense response 
Defense response 
Defense response 

Stress response 



Physiological function 
sub-category 

ROS detoxification & 
degradation 
Pathogenesis-related 
protein 

ROS production & 
hypersensitive response 
induction 
Antioxidative and 
peroxidase activity 



Diff. ab. Accession No. 



53 
52 
33 

84 



GSTUM_ARATH 
MTR_2g035100 
RBOHB_SOLTU 

MTR_4g132110 



153 


3V 


Heat shock protein 83 


Stress response 


Chaperone activity 


73 


HSP83JPONI 


923 


24V 


Peroxidase pxdC 


Stress response 


Antioxidative and 
peroxidase activity 


68 


Q40366_MEDSA 


129 


3V 


Chloroplast envelope 
membrane 70 kDa 
HSP-related 


Stress response 


ATP binding 


45 


HSP7E_SPIOL 


396 


3V 


Stress responsive 
alpha-beta barrel domain 
protein 


Stress response 


Uncharacterized protein 


35 


TC177020 


735 


3V 


Nodule-specific 
cysteine-rich peptide 96 


Signal transduction 


Nodulation 


75 


Q2HW73_MEDTR 


417 


3V 


CBL-interacting protein 
kinase 15 


Signal transduction 


Hormone metabolism 


56 


CIPKF_ORYSJ 


718 


5V 


Auxin responsive SAUR 
protein 


Signal transduction 


Hormone metabolism 


47 


TC 196766 


266 


3V 


Dihydroflavanol-4- 
reductase 1 


Secondary 
metabolite 
biosynthesis 


Flavanoid biosynthesis 


34 


Q6TQT1_MEDTR 


891 


3V 


H+-transporting two-sector 
ATPase, C (AC39) subunit 


Transport 


Ion transport 


49 


MTR_7g009590 


MtROP9i. 


G. INTRARADICES 










6 


5R 


Class-10 

pathogenesis-related 
protein 1 


Defense response 


Pathogenesis-related 
protein 


6 


PR1_MEDSA 


10 


5R 


Pprg2 protein 


Defense response 


Pathogenesis-related 
protein 


6 


Q8L6K8_MEDSA 


52 


3R 


Haem peroxidase, 
plant/fungal/bacterial 


Defense response 


ROS detoxification & 
degradation 


5 


MTR_7g086820 


111 


3R 


MLP-like protein 


Defense response 


Pathogenesis-related 
protein 


4 


MTR_8g045570 


408 


5R 


Endochitinase 


Defense response 


Cell wall degradation 


3 


TC1 73292 


14 


3R 


Peroxidase2 


Stress response 


Antioxidative and 
peroxidase activity 


14 


O24080_MEDSA 


10 


3R 


Heat shock 70 kDa protein, 
mitochondrial 


Stress response 


Chaperone activity 


10 


HSP7M_PEA 


371 


5R 


Prohibitin 1-like protein 


Stress response 


Chaperone activity 


4 


Medtr4g078200.1 


707 


3R 


Clathrin assembly protein: 
At4g32285 


Stress response 


Endocytosis 


4 


CAP1_ARATH 


422 


5R 


lnositol-3-phosphate 
synthase 


Signal transduction 


Protein-protein 
interaction 


23 


IN01_NICPA 


496 


24R 


GTP-binding protein-related 


Signal transduction 


GTP binding / GTPase 
activity 


13 


Q9AY71_0RYSJ 


402 


24R 


Soluble inorganic 
pyrophosphatase 


Signal transduction 


Protein-protein 


12 


Medtr2g01 0430.1 








interaction 
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Table 2 | Continued 



Spot No. 


hpi 


Protein name 


Physiological 
function 


Physiological function 
sub-category 


Diff. ab. 


Accession No. 
















423 


24R 


NAD(P)H-dependent 


Secondary 


Phytoalexin biosynthesis 


24 


6DCS_SOYBN 






6'-deoxychalcone synthase 


metabolite 














biosynthesis 








570 


24R 


Tic22 


Transport 


Intracellular transport 


14 


Medtr3g101 630.1 


Mtvector. 


G. INTRARADICES 










219 


3V 


Cationic peroxidase 1 


Defense response 


ROS detoxification & 


43 


PER1_ARAHY 










degradation 






747 


24V 


Pathogenesis-related 


Defense response 


Pathogenesis-related 


20 


Q9M7D9_PEA 






protein 4A 




protein 






761 


5V 


Glutathione peroxidase 


Defense response 


ROS detoxification & 


23 


Medtr8g105630.1 










degradation 






570 


5V 


Polygalacturonase inhibitor 


Defense response 


Cell wall degradation 


14 


Medtr7g023590.1 






protein 










862 


5V 


Pathogenesis-related 


Defense response 


Pathogenesis-related 


13 


MTR_2g012370 






protein 1a 




protein 






262 


3V 


Peroxidase2 


Stress response 


Antioxidative and 


23 


O24080_MEDSA 










peroxidase activity 






22 


5V 


Thioredoxin hi 


Stress response 


Antioxidative and 


22 


A1BLP6_MEDTR 










peroxidase activity 






41 


5V 


Prn+Gin P?1 l^inaco inhihitnr 
nuicii I I Ml IQoO 11 11 IIUI IKJ\ 


Stress response 


Signal transduction 


17 


P21_SOYBN 










inhibitor 






96 


24V 


Small heat shock protein, 


Stress response 


Heat shock protein 


10 


HS21C_Wheat 






chloroplastic 










772 


5V 


Heat shock protein DnaJ 


Stress response 


Chaperone activity 


4 


MTR_7g114150 


477 


3V 


14-3-3-like protein 


Signal transduction 


Protein-protein 


34 


Medtr3g099380.1 










interaction 






745 


5V 


Caffeic acid 


Signal transduction 


Protein-protein 


14 


Medtr4g038440.1 






3-O-methyltransferase 




II lie! dU UUI 1 






545 


3V 


Blue (Type 1) copper 


Signal transduction 


Ion binding 


12 


MTR_7g086090 






domain 










725 


3V 


Progesterone 


Secondary 


Steroid biosynthesis 


16 


Medtr3g01 3890.1 






5-beta-reductase 


metabolite 














biosynthesis 








112c 


5V 


v-H(+)-ATPase subunit A 


Transport 


Ion transport 


39 


D7EYG6_SOYBN 



List of proteins identified with highest induction in abundance for the selected physiological function categories after infection with S. meliloti, A. euteiches and G. 
intraradices for the time points 1 h, 3h, 5h, and 24 h. Proteins names that are underlayed with gray color are discussed in the manuscript. 



clear visualization of induction patterns of all identified proteins 
according to differences in abundance. The presented Heatmap 
(Figure 2; https://gelmap.de/532) was created by clustering all 
protein induction values of proteins of similar physiological func- 
tions both for major and sub-categories found at a certain selected 
infection time-point. 

Overall, a comparative analysis of induced protein patterns 
considering all infections showed a distinctively lower num- 
ber of total proteins induced in MtROP9i transgenic roots 
both for pathogenic and symbiotic infections as compared to 
Mtvector roots (Figure 2; Supplementary Figures 1A-C, 2A-C). 
This indicates the responsive role of ROP GTPase MtROP9 
protein during early microbial infection signaling. Hereby - 
under conditions of hampered ROS synthesis — induction of 
early defense-related infection protein network mainly com- 
prised of ROS-related enzymes was significantly affected in 



MtROP9i roots, forming the first line of plant host defense 
(Figure 2; Supplementary Figures 1A-C). The amount of total 
protein induction during rhizobial infection was found signif- 
icantly higher as compared to the two fungal infections inves- 
tigated. Considering fungal infections, the induction pattern 
reached maximum at 5 hpi in control roots with significant 
reduction at 24 hpi, but with rhizobia infection, the induction 
reached maximum as early as at 3 hpi remaining constant even 
at a later timepoint (Supplementary Figures 1A-C: minigraphs; 
Supplementary Figures 2A-C). 

Visual evaluation of this Heatmap-GelMap for microbial 
infections in MtROP9i and Mtvector led to detection of major 
results as exemplarily summarized in the following: 

- The total number of induced proteins identified in MtROP9i 
transgenic roots after symbiotic and pathogenic infections were 
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Physiological function Type of Infection: Sinorhizobium meliloti 

main category u. trvncatula line MtROP9i Mtvector control 
Timepoint 'R 3R 5R24R 1V 3V 5V 24V 



defense response 



Aphanomyces euteiches Glomus intraradices 

MtROP9i Mtvector control MtROP9i Mtvector control 

1R 3R 5R 24R 1V 3V 5V 24V 1R 3R 5R 24R 1V 3V 5V 24V 



Physiological function 
sub-category 



t 3K t>K 24K IV JV SV 24V 



stress response 



signal transduction 



energy metabolism 



protein folding & 
processing 



amino acid metabolsim 



nucelic acid metabolism 



secondary metabolite 
biosynthesis 



ion binding & cofactor activity 



transport 



lipid/fatty acid metabolism 



cell development & division 



nutrient substrates storage 



hormone metabolism 



structural protein 
orn.er. development 




The heatmap legend 

Heatmap of total protein induction in transgenic roots of two different Medicago truncatula (Mt) populations at early infection time-points 
with S. meliloti, A. euteiches and G. intraradices. Scores represent the induction fold in one Mt population as compared to the other. 

- MtROP9i: Mt transgenic roots derived from RNAi gene silencing of GTPase MtROP9 {TC173331-MtGI), N = 12 populations of 200 plants 

- Mtvector: Mt. transgenic roots derived from introduction of empty vector construct. N= 12 populations of 200 plants 

- Timepoints: 1 R - 1 hour after infection in MtROP9i roots, 3R, 5R and 24R, respectively 

1V - 1 hour after infection in Mtvector roots, 3V, 5 V and 24V, respectively 

- Protein regulation/induction ranges from: 

Light red colour code - Low (>1 .5fold) induction r — 

Dark red colour code - High (<285 fold) induction j , 



FIGURE 2 | The heatmap legend. 
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significantly lower as compared to control roots, indicating that 
the infection proteome network was affected after silencing the 
signaling protein MtROP9 and therewith the perpetual sup- 
pression of ROS defense signaling. Interestingly, the individual 
protein inductions following rhizobia infection with S. meliloti 
was the highest in control roots as compared to fungal infec- 
tions with pathogenic oomycete A. euteiches and symbiotic 
interactions with G. intramdices. The high protein inductions 
in rhizobia infections could be attributed to the high specificity 
of S. meliloti — M. truncatula infection and genomic adaptation 
of the host plant, but also to the similarities between rhizo- 
bial bacteria and pathogenic bacteria as shown via phylogenetic 
analyses based on their sequences (Willems and Collins, 1993; 
Young et al., 201 1). In addition, several other proteins have been 
induced to promote formation of rhizobia symbiosis such as 
the nodule-specific cysteine-rich peptide 96, dihydroflavanol- 
4-reductase 1. However, the total number of individual proteins 
induced in A. euteiches infections was higher, where it can be 
surmised that the host in addition to induction of the early 
ROS defense signaling enzymes engages other alternative path- 
ways particularly at the advanced stages of infection process, 
including induction of cell wall degrading enzymes such as 
polygalacturonase inhibitor, pectinesterase and SGT1 protein 
as well as PR proteins. 

- Induction of proteins involved in ROS production and 
detoxification was significantly reduced in MtROP9i roots 
as compared to control roots after microbial infection 
(Supplementary Figures 2A-C). RBOH 1 and RBOH B were 
highly induced at the early points of S. meliloti infection par- 
ticularly in control roots but clearly reduced in MtROP9i. In 
addition, their induction levels considering G. intramdices and 
A. euteiches infections were moderate in control roots. SOD 
[Cu-Zn] was highly induced after S. meliloti infection, but its 
induction was noted only for the three timepoints 3, 5, and 
24hpi. The decline in induction of proteins responsible for 
ROS production and the concomitant decrease of detoxifying 
enzymes at 12 hpi is indicative of the essential role of ROS dur- 
ing root and nodule development (Lohar et al., 2007; Puppo 
et al., 2013). A transient decline in ROS activity is thought to 
be essential for early NF signaling and the infection process at 
the root hairs (Lohar et al., 2007). In contrast, permanently 
suppressed levels of ROS due to inactive RBOH proteins as 
documented for MtROP9i transgenic roots led to a drastic 
impairment of early rhizobial infection accompanied by abnor- 
mal root hair deformations (Kiirika et al., 2012). Regained 
accumulation of ROS after 12 hpi in control roots could have 
multiple functions, one as part of a typical defense response 
to limit bacteria entry, second, as compounds needed for pro- 
gression of infection thread or as signals for symbiotic protein 
synthesis (Soto et al., 2011). 

- 92 proteins, classified as stress-related proteins were highly 
induced after S. meliloti infection with 34 and 58 proteins 
identified in gels of MtROP9i and control roots, respec- 
tively. Among them, the peroxidases (84 fold at 5 h) involved 
in scavenging of peroxide and detoxification of ROS was 
highly induced. Several other enzymes including peroxidases, 
catalases, SOD (Cu-Zn) as well as enzymes involved in ROS 



production and hypersensitive response induction such as 
harpin binding protein 1, hydroxyacylxylotathione hydrolase 
as well as perforin domain containing protein 1 were also 
detected. 

- 109 proteins involved in signal transduction were found to 
be highly induced where 33 and 76 proteins were induced 
in MtROP9i and controls, respectively. Calmodulin (CaM) 1, 
2, 8 and CaM binding proteins were induced in the cate- 
gory of calcium second messengers, which constitutes pro- 
teins involved in the initial plant-microbe interactions. CaM 
mediates the calcium-dependent signaling by functioning as 
a decoder for the Ca 2+ signatures during signal transduction 
especially at the early points of infection (Bender and Snedden, 
2013). Induction of CaM depicts progressive colonization of 
the tissue by the pathogen thereby eliciting the secondary line 
of defense barrier. Calreticulin protein that functions in con- 
trolling plant defense by regulating the concentration of Ca 2+ 
ions in the cell during signaling pathways as well as acting 
as molecular chaperone (Qiu et al, 2012) was induced only 
in control roots, particularly in fungal infections and not in 
MtROP9i. It was recently shown to confer resistance against 
oomycete pathogen phytopthora infestans in Nicotiana ben- 
thamiana (Matsukawa et al, 2013). Interestingly, the CaM and 
calreticulin proteins were highly induced in the Mtvector roots 
which also had concomitant high induction of ROS-related 
enzymes, suggesting a connection between ROS and calcium 
signatures. Both signals most likely represent a cross-talk that 
is constituted to modulate downstream nuclear activity result- 
ing in induction of pathogenic or symbiotic specific protein 
networks. 

- Proteins involved in antioxidative and peroxidase activity such 
as the thioredoxins were induced mainly at 3 and 5 hpi, espe- 
cially after A. euteiches infection. The cytosolic form of thiore- 
doxin h functions in response to oxidative stress which may 
have occurred at the early points of microbial infection, mainly 
pronounced in the control roots. Thioredoxin also accumulates 
in self-incompatibility reactions, seed germination and early 
seedling development (Chi et al., 2013). 

- Paltry induction of ROS-related enzymes in MtROP9i trans- 
genic roots after infection could not be ruled out especially 
after S. meliloti infection. This could be attributed to a knock- 
down but not a fully knockout of MtROP9 proteins or even 
the presence of other related GTP binding protein partners that 
as a consequence could potentially contribute to ROS synthesis 
at the cellular plasma membrane, even though at significantly 
reduced levels as compared to control roots leading to overall 
impairment of rhizobial infection signaling. Alternative pro- 
teins present in MtROP9i which belong to the category GTP 
binding/GTPase activity include the dynamin-related protein 
and guanine nucleotide -binding proteins. 

- Proteins specific for rhizobia infection, such as the nodule- 
specific cysteine-rich peptide responsible for early nodulation 
signaling, were highly induced at 3 hpi after S. meliloti inocu- 
lation exclusively in control roots. The predominant induction 
of this protein pattern very well indicates the onset of rhizo- 
bial infection signaling in control roots. Since proteins involved 
in early infection signaling were not detected in MtROP9i 
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roots, evidence for a hampered rhizobial infection is given as 
previously reported (Kiirika et al., 2012) 

- G proteins have been shown to play a key role in early nodu- 
lation signaling (Choudhury and Pandey, 2013). Hence, low 
induction of proteins involved in nodulation in MtROP9i pro- 
vides evidence to the interfered signal transduction process 
during root-rhizobia symbiosis due to MtROP9 GTPase silenc- 
ing as one of the major GTPases for early nodulation signaling. 

- Major latex protein (MLP) and ABA-response protein (ABR17) 
related to the PR- 10 proteins category were predominantly 
induced after A. euteiches and G. intmrctdices infections at 
5 hpi in control as compared to MtROP9i. These proteins are 
induced when the intracellular ABA concentration increases 
in the cell during the course of infection. Generally, PR 
proteins comprising class 10, 1A, 5-1, and 1 as well as 
osmotin/thaumatin-like proteins were induced particularly in 
control roots infected with S. meliloti, indicating an introduc- 
tion of a profound defense response also in rhizobial infections 
when ROS as signaling molecule is formed. 

- Induction of PR proteins was previously reported both at the 
transcript and protein level (Colditz et al, 2004, 2005, 2007) 
forming the major components of molecular host defense 
response on M. truncatula after A. euteiches infection both at 
the early and later infection stages. The induction of PR pro- 
teins, in particular, was detected distinctively at advanced stages 
of A. euteiches infection (5 and 24 h) with predominantly high 
inductions in control roots as compared to MtROP9i roots. So 
far, reports indicate that 14 classes of PR protein are known 
in plants (Spoel and Dong, 2012). Previously, it was reported 
that especially both PR protein classes 5 and 10 were con- 
jointly modified in expression following A. euteiches infections 
in M. truncatula (Colditz et al, 2007; Trapphoff et al, 2009). 
Several other proteins differentially induced were classified in 
the category of defense response, proteins involved in pro- 
tein degradation, proteolysis and proteolysis inhibitor, induced 
only after S. meliloti and A. euteiches. Proteins involved in pro- 
teolysis, chaperone activity, protease inhibitor (kunitz-type) 
and protein modification were found induced exclusively after 
A. euteiches infection in control roots. 

- In the category of secondary metabolite and biosynthesis, 
proteins responsible for phenylpropanoid biosynthesis 
(e.g., hydroxycinnamoyl-CoA quinate hydroxycinnamoyl- 
transferase), phytoalexin biosynthesis (e.g., chalcone-fiavonone 
isomerase 1, isoflavonoid xylocosyltransferase and NAD(P)H- 
dependent 6'-deoxychalcone synthase) and flavonoid 
biosynthesis (e.g., dihydroflavanol-4-reductase 1 and CXE 
carboxylesterase) were found induced mainly in control roots. 

- Structural protein (e.g., MFP1 attachment factor), hormone 
metabolism (e.g., abscisic acid receptor PYR1 and cytokinin-O- 
glucosyltransferase), cell division and development (e.g., actin, 
profiling and cyclin) and lipid/fatty acid metabolism (e.g., 
phospholipase D, fatty acid oxidation complex subunit alpha 
and lipoxygenase) were highly induced in control roots espe- 
cially in Rhizobia infection as compared to fungal infections. 

- Proteins playing key role in nodulation, sulfur metabolism 
and pyrimidine biosynthesis were only induced with S. meliloti 
infection in control roots. 



- Proteins responsible for protein-protein interaction (14-3-3 
like protein), ion binding (blue type copper domain) and 
NAD binding were only induced after G. intraradices infections 
(Table 1). 

DEFENSE PROTEOME OF M. TRUNCATULA MtR0P9l TRANSGENIC 
ROOTS DEFECTIVE IN ROS SIGNALING 

Of the 17 main physiological functional categories of induced 
proteins identified during the evaluations, significantly high 
numbers of proteins were found to be involved in defense 
response, stress response, signal transduction and energy 
metabolism. 

During S. meliloti infections, 92 defense-related proteins were 
found to be induced with only low protein number (27 proteins) 
induced in MtROP9i as compared to control roots, where 65 
proteins were induced. The majority of defense-related proteins 
induced in MtROP9i roots were those involved in alternative 
defense response pathways and not linked to ROS, indicating that 
the potential for host cell to synthesize ROS was compromised. 
Induction of protein SGT1 homolog (77 fold at 3 h) involved in 
protein degradation and the pathogenesis-related protein PR10-1 
was found to be significantly high in MtROP9i roots as com- 
pared with Mtvector roots. Mitochondrial prohibitin (12 fold at 
3 h) with the role of maintaining the integrity of the organelles 
was also significantly induced in MtROP9i. Recent findings shows 
that prohibitins are involved in mediating stress tolerance (abiotic 
stress, pathogen infection and elicitor signaling) as well as trigger- 
ing retrograde signals in response to mitochondrial dysfunction 
(AkenetaL, 2010). 

Signaling proteins, such as proteins of the germin-like protein 
subfamily 1 and 2 (41 fold at 5h and 45 fold at 24 h), which are 
carbohydrate binding proteins, were highly induced in MtROP9i 
roots. This indicates that early nodulation signaling is gener- 
ally not affected by the gene silencing. Since the nodule-specific 
cysteine-rich peptide 96 (78 fold at 5 hpi) which plays a key role 
in the early infection signaling, was highly induced in control 
roots but not found in MtROP9i roots, it is much likely that infec- 
tion processes shortly after the NF signaling are ROS-dependent 
via cross-linking of attachment proteins. Thus, in case of ROS 
depletion, lacking attachment of bacteria to root hairs results in 
impairment of rhizobial infection processes. 

During A. euteiches infections, 165 proteins involved in defense 
response were identified, of which a significantly low number 
of proteins (64 proteins) were found induced in MtROP9i as 
compared to controls (101 proteins). Majority of these proteins 
induced in MtROP9i roots comprised those involved in alter- 
native or secondary defense pathways giving an indication that 
the RBOH activity in ROS biosynthesis was affected. Cell wall 
degradation proteins including endochitinases (33 fold at 5 h) and 
polygalacturonase inhibitor protein (22 fold at 24 h), proteoly- 
sis such as nectarin IV (10 fold at 3 h), PR 5 protein (12 fold at 
3 h) were highly induced in MtROP9i. PR-10-type proteins anno- 
tated as abscisic acid responsive proteins (AB17s) (7 fold at 5h) 
were also induced in MtROP9i roots. Abscisic acid in Arabidopsis 
guard cells was shown to enhance levels of ROS (Pei et al., 2000). 
The meager traces of ROS-related enzymes detected at later stages 
of infection may have been enhanced by the presence of abscisic 
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acid in the infected cells whereof indirectly affirmed by the sig- 
nificant induction of abscisic acid responsive proteins. Induction 
of PR proteins occurred at advanced stages of infection represent- 
ing a secondary line of host defense system especially in MtROP9i 
roots. Their induction was also reported earlier with A. eute- 
iches infections as the major component in the pathogen defense 
response of M. truncatula both at transcript and protein levels 
(Colditz et al, 2004, 2005, 2007; Schenkluhn et al, 2010). 

Interestingly, a cytoplasmic redox protein non-expressor of 
PR1 (NPR1) was suppressed in MtROP9i roots and induced 
only in control roots at 5 hpi after pathogenic infection with 
A. euteiches, implying that its induction was hampered with the 
knockdown of MtROP9 protein. However, it was not induced 
in any of the symbiotic interactions. NPR1 protein functions in 
sensing salicylic acid (SA) during microbial infection where it 
exists in the cytoplasm as an oligomer in its non-induced forms 
but reduced to its monomeric form after pathogen infection and 
translocated to the nucleus leading to induction of PR proteins 
(Rochon et al, 2006; Vlot et al, 2009). Plants produce high con- 
centrations of SA during systematic acquired resistance (SAR) 
both in infected and non-infected tissues that functions in the 
signaling of defense responses including induction of PR proteins 
(Oide et al., 2013). Presence of cytoplasmic ROS was shown to 
inhibit the action of NPR1 protein hence suppressing induction of 
PR proteins in M. truncatula seedlings infected with an incompat- 
ible bacteria Pseudomonas syringae (Peleg-Grossman et al., 2012). 
Based on our analyses, the suppressed induction of NPR1 pro- 
tein at the early points of infections in MtROP9i but its particular 
induction in control roots could be due to high levels of ROS 
production in Mtvector. 

The predominant induction of glutathione-S-transferase 
(GST) in control roots infected with A. euteiches with minimal 
induction in MtROP9i roots indicates the reduced activity of GST 
redox system responsible for primary antioxidative processes in 
the latter. The knockdown of MtROP9 also affected the expression 
of stress-related proteins. 51 proteins were induced in MtROP9i 
which was less compared with control roots where 81 proteins 
were induced after A. euteiches infection. Overall, 132 stress- 
related proteins were identified, which exhibited high inductions. 
Majority of stress-related proteins induced in MtROP9i roots 
were not linked to ROS production as found in control roots, 
but form part of other alternative cellular reaction pathways. The 
protein P21 kinase inhibitor, a potent cyclin- dependent kinase 
inhibitor (CDKI) was highly induced (50 fold at 3 h). P21 kinase 
inhibitor is induced as stress reaction of plants by deregulation of 
cell proliferation controlling growth pathways cell-cycle regula- 
tory and integration of developmental signals in the cell machin- 
ery (De Veylder et al, 2001). A glycine-rich RNA-binding (3 fold 
at 24 h) as well as the cold shock protein- 1(46 fold at 3h) were 
induced in MtROP9i and are involved in post-transcriptional reg- 
ulation of gene expression in plants under biotic or abiotic stresses 
(Mangeon et al, 2010). 

Within G. intraradices infections, 43 defense-related proteins 
were induced where 15 and 28 proteins were induced in MtROP9i 
and control roots, respectively. PR proteins of class 10 (78 fold at 
24 h) and Pprg2 (6 fold at 5 h) were induced in MtROP9i roots. 
The major latex protein (MLP) was detected both at 3 and 5 hpi. 



MLPs play a role in host defense (Betsy et al., 2009) and were only 
detected in fungal infections potraying their expression speci- 
ficity to fungal infections. 92 stress-related proteins were highly 
induced after G. intraradices infection where 42 and 50 proteins 
were found in MtROP9i and control roots, respectively. Among 
the induced proteins in this category, clathrin assembly protein 
(At4g32285) (4 fold at 3 h), stress related proteins involved in the 
formation of vesicles in the cytoplasm for intracellular protein 
trafficking was highly induced in MtROP9i. Majority of highly 
induced proteins detected after G. intraradices infection in control 
roots comprised those specific for this type of infection, including 
thel4-3-3-like protein (34 fold at 3h), a family of signaling pro- 
teins, blue (type 1) copper protein (29 fold 5 h) involved in signal 
transduction as well as v-H(+)-ATPase subunit A (39 fold 5 h) 
playing a key role in membrane transport. 69 proteins involved in 
signal transduction were induced with 23 and 46 proteins induced 
in MtROP9i and control roots, respectively. The soluble inor- 
ganic pyrophosphate (12 fold at 24 h) involved in protein-protein 
interactions during signal transduction was highly induced in 
MtROP9i as compared with controls. 

Several other proteins were differentially induced such 
as translationally-controlled tumor protein homolog (TCTP) 
exhibiting low induction in MtROP9i especially after A. eute- 
iches infection as compared to control roots, which suggests the 
influence of MtROP9 knockdown. TCTP is shown to be a Ca 2+ 
-binding protein ubiquitously expressed in all eukaryotic cells 
(Zhang et al, 2013). We reported its high induction at the tran- 
script level in tomato treated with elicitors silicon and chitosan 
against soilborne bacteria Ralstonia solanacearum (Kiirika et al., 
2013). Recent studies have also shown its involvement in negative 
regulation of hypersensitive reaction in N. benthamiana plants 
challenged with R. solanacearum (Gupta et al., 2013). Two modes 
of action of TCTP cytoprotective activity have been suggested, 
where it is said to act as a sequester of Ca 2+ hampering pro- 
grammed cell death (PCD) by reducing levels of Ca 2+ in the 
cytosol or interact with other cytosolic membranous proteins of 
the cytosolic PCD machinery, thereby mitigating the progress in 
cell death. Caffeic acid 3-O-methyltransferase (CCOMT), a pro- 
tein involved in conversion of caffeoyl-CoA to synapoyl-CoA, 
an intermediate in lignification in cells under stressed condition 
(Kosova et al, 2013) was induced particularly with the fun- 
gal infections while exhibiting low inductions in MtROP9i as 
compared to Mtvector roots. Taken together, the mechanisms of 
fungal invasion and ROS scavenging to the cell triggers CCOMT 
synthesis to function putatively as a compensatory response by 
enhancing lignin formation so as to reinforce the cell wall and 
maintain its functional integrity. 

CONCLUSIONS 

The early molecular events in the cell upon microbial infection 
comprises ROS defense signaling as major first defense barrier. 
ROS generation is mainly related to activation of RBOH pro- 
teins which is achieved by protein interaction with cellular small 
GTPases, thus functioning as molecular switches. Therefore, from 
the current investigations, we have shown that silencing of small 
GTPase MtROP9 in M. truncatula roots results in significant sup- 
pression of ROS-related enzymes at the early points of infection 
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both for symbiotic and pathogenic interactions and a concomi- 
tant activation of alternative defense pathways by the host cell to 
counteract the infection progress. The induction of ROS-related 
enzymes such as RBOH 1/B, glutathione-S-transferase, peroxi- 
dases, catalases, SOD was found clearly suppressed in MtROP9i 
transgenic roots as compared to control roots. This implies 
that the RBOH activity responsible for oxidative burst during 
microbial infection was affected. Instead, proteins involved in 
alternative defense pathways such as cell wall degradation, chitin 
hydrolysis, proteolysis, and PR protein induction were highly 
induced. Other proteins whose induction was affected also were 
involved mainly in signal transduction pathways, general stress 
response, and chaperone activities. Since proteins involved in 
energy metabolism exhibited minimal differential induction in 
MtROP9i as compared with control roots, it clearly shows that 
ROS signaling is the major force for introduction of effective 
defense signaling, while enzymes related to energy metabolism 
which represents a prerequisite for a profound host infection 
defense generally remain retrievable. The total amount of pro- 
tein induction was highest after S. meliloti infection in Mtvector, 
followed by the induction after pathogenic A. euteiches infection. 
Despite the overlapping similarities existing between pathogenic 
and symbiotic interactions, the host plant has to maintain speci- 
ficity such that it activates response cues that are either meant 
to promote or desist the interaction. This process consists of 
broad array of protein networks forming central dynamics in the 
cell, which when analyzed give huge proteome datasets. The new 
Heatmap-GelMap module presented allowed systematic and easy 
visualization of all induced proteins displaying also their physi- 
ological functions in the cell. To this end, the elaborate protein 
induction evaluation upon symbiotic and pathogenic interac- 
tions reported herein, depicts the essential role of small GTPase 
MtROP9 during the early ROS defense signaling in plants. 
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SUPPLEMENTARY MATERIAL 

The Supplementary Material for this article can be found 
online at: http://www.frontiersin.org/journal/10.3389/fpls.2014. 
00341/abstract 

Supplementary Figure 1 | Protein induction patterns of M. truncatula 
MtROP9i and Mtvector roots after S. meliloti (A), A. euteiches (B) and G. 
intraradices (C) infection visualized on a Heatmap. The heatmap shows 
distinct proteome subsets created from total induction values of all 
identified proteins classified into different physiological functions. The 
color intensity potrayed as a shift in gradient from light red to deep red 
represents the relative changes from a minimum to a maximum value 
allowing easy spotting of trends or standouts in each infection at a given 
time after infection. The total induction for each category is shown on the 



right and the minigraphs for the trend of total induction over 1 , 3, 5, and 
24 hpi is shown at the bottom. 

Supplementary Figure 2 | Heatmap of total protein regulation in each 
physiological category over the four time-points for S. meliloti (A), 
A. euteiches (B) and G. intraradices (C) infections. 
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